Stress-energy of a quantized scalar field in static wormhole spacetimes 
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An analytical approximation of {Tfj,^) for a quantized scalar field in a static spherically symmetric 
spacetime with a topology x is obtained. The gravitational background is assumed slowly 
varying. The scalar field is assumed to be both massive and massless, with an arbitrary coupling ^ 
to the scalar curvature and in a zero temperature vacuum state. It is demonstrated that for some 
values of curvature coupling the stress-energy has the properties needed to support the wormhole 
geometry. 
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In recent years Lorentzian wormholes have been of some interest as nontrivial geometric and physical objects. 
The topology of such spacetimes is x S^. The characteristic feature of wormhole geometry is the existence of 
a sphere of minimal square which cannot be contracted to the point. The possibility of the existence of a static 
spherically symmetrical traversable wormhole as a topology-nontrivial solution of the Einstein equations was first 
' studied by Morris and Thorne [Q. They found that the matter which threads the wormhole throat should have 
' unusual properties. In particular, the radial tension of the matter should exceed its density both locally at the throat 
101 and integrally along the radial direction ||^. As an example of such matter they suggested the Casimir vacuum 
^ [ between conducting spherical plates. Violations of the energy conditions in static wormholes have been analyzed in 
\^ • detail by a number of authors (see, for example, Q). Recently Barcelo and Visser demonstrated that a scalar 
\^ ' field with a positive curvature coupling violates all the standard energy conditions. They found the entire branch 
T-H of traversable wormhole solutions for gravity plus a nonminimally coupled massless scalar field. Solutions describing 
On ' primordial wormholes are obtained for A/" = 4 SU(N) super Yang-Mills theory in [^. The possibility of the existence of 
I wormhole solutions has been investigated in Brans-Dicke theory [Q-^ , higher derivative gravity , and Gauss-Bonnet 
g ' theory 0. 

Within the framework of general relativity as an example of the matter that supports the wormhole geometry one 
. can consider the vacuum of quantized fields. This approach gives the possibility of describing the wormhole metric as a 
"V* ' self-consistent solution of the semiclassical Einstein field equations ||^,|l3| . The principal problem encountered in this 
Oh, method is the impossibility of calculating the functional dependence of the stress-energy tensor of quantized fields on 
^ • the metric tensor. Nevertheless, in some cases approximate expressions can be obtained [p^ -p9[. In above mentioned 
work |l^,|l^ the approximations of Frolov-Zel'nikov and Anderson-Hiscock-Samuel [ p9[ were used. But as noted 
in pT[ | in spacetime that is a direct product of the Minkowski plane and a two-dimensional sphere of fixed radius these 
approximations are not applicable. The approximate expressions for the electromagnetic and neutrino stress-energy 
tensor obtained by Khatsymovsky |^,^ do not have this deficiency. In this work the results were obtained by the 
WKB method, and, as the zero order of the expansions was used, the exact solutions of the mode equations were 
obtained in spacetime with the metric ds^ = —fadt^ + dP -f rQ{d9^ -t- so? 9 d^p^) (/o and tq are constants), i.e., in 
the above mentioned spacetime. It is well known that the vacuum expectation value of the stress-energy tensor of 
quantized fields depends on the topology of spacetime. In this connection the approach of Khatsymovsky seems more 
natural in wormhole spacetimes. 

In this paper an approximate expression for the stress-energy tensor of a quantized scalar field in static spherically 
symmetric spacetimes with topology E? x is calculated. The Anderson-Hiscock-Samuel approach 129] is used but 
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the zeroth order of the WKB solution of the mode equation is redefined by analogy with 1 26 2^] and mode sums are 
computed exactly as in pO| , ^ . 

The units Ti = c = 1 are used throughout the paper. 
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II. AN UNRENORMALIZED EXPRESSION FOR {Tj^) 



In this section the Euchdean space approach is used to calculate an unrcnormalized expression for {Tj^) for a scalar 
field in a static spherically symmetric spacetime with topology 5^ x . The metric for a general static spherically 
symmetric spacetime when continued analytically into Euclidean space can be written as 



ds^ = f{p)dT^ + dp^ + r^{p){de^ + sin^ddip^). 



(1) 



We assume that — cxi < p < oo, / and r are arbitrary functions of p, and r is the Euclidean time (r = it, where t is 
the coordinate corresponding to the timelike Killing vector which always exists in a static spacetime). 

(r^) of a quantized scalar field (f) can be computed using the method of point splitting from the Euclidean Green's 
fimction Ge{x,x) as follows 



(T,^)™re„ = (1/2 - (5''"G£;fi, + 9:Ge;^s.) + (2^ - l/2)5^:.g^"GB;,fi - ^e.^u 



(2) 



where m is the mass of the scalar field, ^ is its coupling to the scalar curvature i?, and g'^ is the bivector of parallel 
transport of a vector at x to one at x. 

In 1^ the form of Ge{x, x) was derived for a scalar field in a static spherically symmetric spacetime when the field 
is in the zero temperature vacuum state defined with respect to the timelike Killing vector 



GEix;x) = / dtJCOs[w(T - f)] ^(2Z + l)P;(cos7) C^/ p„/(/9<) g^;(/9>), 

^ 1=0 



(3) 



where Pi is a Legendre polynomial, cos 7 = cos 6* cos 6* + sin sin 6* cos (93 — <f), C^ji is a normalization constant, /3< and 
/9> represent the lesser and greater of p and p, respectively, and the modes p^^i (p) and q^^i (p) obey the equation 



dp^ 



1 df 1 dr^ 



2/ dp dp 
They also satisfy the Wronskian condition 



dp 



dq^ 
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= 0. 
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(4) 
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After the substitution 



1 



1 



exp 



: exp 



\/2r^W 

it is easy to see that the Wronskian condition (|^) is obeyed if 

C^i = 1 

and the mode equation gives the following equation for W{p) 



Wf-^'^dp 



Wr^/'^dp 



(6) 



W'^^uj'^ + l^ + 1) + 2^ + mV^] 



fiw^y f (wy 5f{w^y 
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V, 



where 
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4^4 



r 

f 



2\n 



2/2 



(r2)'2 /'(r2)' 



(7) 



(8) 



(9) 



The prime denotes the derivative with respect to p. 

Equation (0) can be solved iteratively when the metric functions f{p) and r'^{p) are varying slowly, that is. 
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^WKB = L^,/ L ^ 1, 



(10) 



where 



2f 



-1/2 



and L is a characteristic scale of variation of the metric functions: 



L ^ = max 



1/2 



[Hfr')]' , [Hfr^)]" , [Hfr')] 



1/3 



The zeroth-order WKB solution of Eq. (|8|) corresponds to neglecting terms with derivatives in this equation 
where 



^ / 2 



(11) 



(12) 



(13) 



(14) 



Let us stress that Uq is the exact solution of Eq. (||) in a spacetime with metric ds^ = f odr'^ + dp^ +rQ{d9'^ +sm^ 9dip^), 
where /o and vq are constants. 

Below as in Ref. I^lll it is assumed that 



or 



The second-order solution is 



W=Uo + V- 



> 
2 ^ 1-8^ 

fW fUo" 5fUo'^ 
8 Uo A Uo 16 U§ ' 



(15) 
(16) 

(17) 



Now we can rewrite expression (g) using expressions (||) and (||) and then suppose p = 5^ ^ O.cj) — (j). The superficial 
divergences in the sums over I that appear in this case can be removed as in Refs. [P3,E9|: 



\ t / 1 



2/ 2 ' ■' \ 2f r2 



(9 



tpN2 



B4 



e/ ( -7^^ + rri' + ] {g'n' - ( 2e - ) + ( 2e - ^ ) (m- + ^i?) 



1 \ 1 



4j,2 



/ d 



2 J 4r2 
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(18) 



P I unren 
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2/ r2 



eU7 + ^ (l + (.9'"^)') + ^(.9''*)' 
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^(i + (5'"^n(7^--^-^^)-(2e-7)7^ + (2e-^)(-^ + €i?) 



1 



47- 



1 \ 1 



B,+^^'Jg^fg'\^{^-^B,+^^ 



f d 
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2 y 4r2 
2^*^^ -f 2gPPgP'^ 



1 



(19) 
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where 



(Te) 



0, ^\ "-^ ^ R 



2e 



1 



Si + i2 ( 2e - i ) .9 



j,2 



(20) 



1 /"^ °° 1 

•^0 i=o L 



/ (^ + 1/2) 



(21) 



B. = 



2 /"OO CXD ^ 

— dwcos(ue)^-j 
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2/ 4r2 16/ 
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2 / 16 V r2 
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" / 2 



(22) 



^3 = ^^ '^^cos(ue)^^ 
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16/2 32/ ^2/ 
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(24) 



r(T-r). 



(25) 



III. USE OF THE WKB APPROXIMATION IN EVALUATING (Tif) 



The approximate (of second WKB order) expressions for the quantities Bi, B2, B^, and B4 are obtained by 
substituting the WKB expansion of [Eq. (|)] into Eqs. iM) 



B^ = 



-^^5o(e,M)-^5°(e,M)-Y^ 
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(26) 
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(27) 



B. = 
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+ 



5r4 // 
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16/2 
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(28) 
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The mode sums and integrals in these expressions are of the form 



(29) 



S™{s,fi)= I du cos{eu) 

^0 , n 



a + 1/2) 



2m+l 



;=o [ [w2 + ^2 + (/ + i/2y 



(2«+l)/2 



— subtraction terms 



(30) 



where m and n arc integers, m > and n > —1. The subtraction terms for the sum over I can be obtained by 
expanding the function that is summed in inverse powers of I and truncating at 0(^°)- Such subtracting corresponds 
to removing the superficial divergences in the sums over I discussed above: 



f°° 00 J 

Sq (e, n) = I du cos(eu) < = 
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(Z + 1/2)3 
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+ 



(31) 



^(e,M) = ^"rfwCOs(eu)^|(^« + ^^ ^w2+^2 + (; + i/2) 



(32) 



S'!^{s,n)= I du cos{eu) < 



a + 1/2) 



2m+l 



_ [u2 + ^2 + + 1/2)2 



(2m+l)/2 



(33) 
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For other quantities S^{e, jj) there are no subtraction terms. The details of calculations of S'!^{e, /i) in the limit e ^ 
are discussed in Appendix A: 
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where 
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d 



2 \fj,dfi 



(2to+ 1)!! a 
(2n+ 1) 



^(^-— j 5;^+i(£,Ai) (m>0, n>m + 2), 



°° a;2"-iln|l-a;^| 
1 + e^'^A'^ 



(39) 
(40) 

(41) 



Substitution of these expressions into Eqs. (|2l|)-(|2J) and then into Eqs. ([18D-(|20D gives nontrivial components of 

unren ' 

The rcnormalization of (T^) is achieved by subtracting from (Tl^) unren ^^'^ counterterms (Tl^) jjg and then letting 
f ^ t: 



(T^) = lim [(T^) 

\ u I ren _ u I z 



(42) 



The expression for (^1^) ug are displayed in Ref. [ p^ l and are determined by cr, one half the square of the distance 
between the points x and x along the shortest geodesic connecting them, and cr^, covariant derivative of a. For the 
metric (|h the calculations of these quantities and g^^^ are similar to those in Ref. [M : 
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a' = a* = 0, 
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(47) 



where C is Euler's constant and tods is equal to the mass to ol the field for a massive scalar field. For a massless 
scalar field it is an arbitrary parameter due to the infrared cutoff in {Tj^) j^g. A particular choice of the value of tods 
orresponds to a finite renormalization of the coefficients of terms in the gravitational Lagrangian and must be fixed 
by experiment or observation. 

The procedure described above gives a second-order WKB approximation of (T^^) . The zeroth-order expressions 
of nontrivial components of (T'^)rg„ have the form 
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The second-order expressions of {T^^)^^^^-^ are given in Appendix B. 

Let us stress that the renormalized expectation value of the stress-energy tensor (^8|) , 
const and r{p) = tq = const, i.e., in the spacetime with metric 

ds^ = -Mt^ + V + rlide'^ + sin^ 6*^(^2). 

The expressions (Ha) and (E9h may be simplified in the cases ^ , rn?r'^ ^ 1 and ^ r^, rn^r^ '3> 1. 
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A. The case Lr , m^r^ < 1 



In this case 
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In particular, if ^ = 1/6, 
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In this case tods = "z. Therefore 



B. The case > r^, rn^r^ > 1 
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As was demonstrated by Morris and Thorne ||l|| the radial tension r at the wormhole throat must be positive. In the 
considered case t = — (Tp) > for 

\ P / ren 



(9800 + 735^205) 
210 



1/3 



7 



6(9800 + 735 Vm) 1/3 6 



I w 0.2538. 



(59) 



As an example of a metric function r'^{p) satisfying the condition ^ r^, m^r'^ 1 one can consider a function 



growing as a logarithm: 



ds'^ = -dt^ + dp^ 



EL 

Pa 



(60) 



where ro,a, and po are constants. The conditions ^ and r^m^ ^ 1 gives a^r^/pQ <C 1 and r^m^ ^ 1, 
respectively. 
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IV. CONCLUSION 



We have obtained an analytical approximation of the stress-energy tensor of quantized scalar fields in static spher- 
ically symmetric spacetimes with topology x R^. For some values of coupling to the scalar curvature ^ the 
stress-energy tensor obtained has the needed exotic (in the sense of Morris and Thorne ) properties to support a 
static wormhole. All three approximations (^), (p9|), ( p2[ ) and ( |5^ ) for {Tj^)^^^ are conserved, i. e., 

(-^1/ )ren -.j^i ^ 

and for the conformally invariant field the approximation for {Tl^)^^^ '^^^ ^ trace equal to the trace anomaly. These 
approximations are determined by the correlation between three length scales: L (|l2|), r (|l|) and the Compton length 
1/m. The use of the expression obtained for {Tlf)^^^ as a source term in the semiclassical Einstein field equations 
= SttG {Tl!')ren dcmauds some accuracy. First of all in these equations a new length scale ^/G appears. If on the 
left-hand side of these equations we neglect the terms of order the system of equations will be algebraically 

incompatible. If we keep such terms the terms of corresponding order should be kept on the right-hand side of the 
equations. 

In conclusion let us note that in the case LP' <r^ the WKB parameter Xwkb [Eq. (p^] coinsides with the small 
parameter l/{mL) of the DeWitt-Schwinger expansion of {Tjf) and with the small parameter used in to obtain 
an analytical approximation for {Tjf)^^^. In this case as the first and second WKB orders of (T'^)ren vanish [see Eqs. 



(B1)-(|B3)] and the next orders are needed for calculation of these quantities. 
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APPENDIX A 

The sums over I in expression ( pO| ) can be evaluated by using the following method: 

n n „ 

J„ 2i^/(/ + 1/2) = -27ri^ res [f{z)taTi{Trz)] = - f{z)taTi{Trz)dz, (Al) 
1=0 1=0 ^-'+1/^ J 

where C is a closed contour which surrounds a region on the complex plane containing the poles of tan(7rz): z = ^-|- 1/2 
for I = 0,1, 2, n . It is supposed that countour C is passed anticlockwise, f(z) is a holomorphic function inside this 
contour, and 

\f{x + iy)\ < e{x)e''^y^, < a < 27r, (A2) 

where 

e{x) -^0 for +oo. (A3) 

Let us choose the contour so that 
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where —1/2 < q < 1/2 and h > 0. However, 

pq-\-n-\-l — ih 
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/(z)[tan(7rz + i)]dz 
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where —1/2 < g < 1/2, /(z) is a hofomorphic function for i?ez > q, and /(z) satisfies condition (|A2|). 
Using the last expression we can calculate the sums in Eq. (pG) : 
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Now we can calculate the integrals over u in (|3^) . Note that these expressions must be expanded in 
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3(eu)(M^ + 
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where iir„(x) is Macdonald's function and C is Euler's constant. 
The second type of integral over u has the form 
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If we also take in consideration 
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l + e2'^^ 48' 7o l + e^^^^^ 1920' 
the resulting expresions for Sl^{e,^) can be presented as Eqs. (^-(^. 
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